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IMPINGEMENT OF WATER DROPLETS ON NACA 65,-208 AND 
65--212 AIRFOILS AT 4° ANGLE OF ATTACK 


By Rinaldo J. Brun, Helen M. Gallagher, and Dorothea E. Vogt 


SUMMARY 


The trajectories of droplets in the air flowing past an NACA 654-208 
airfoil and an NACA 65,-212 airfoil, both at an angle of attack of 49, 


were determined. The amount of water in droplet form impinging on the 
airfoils, the area of droplet impingement, and the rate of droplet 
impingement per unit area on the airfoil surface affected were calculated 
from the trajectories and are presented herein to cover the following 
range of conditions: 


Droplet diameter, microns 


Critical flight speed 


Airplane speed, mph 


35 ,000 


Altitude, ft 


Chord. length, ft 20 


The amount, extent, and rate of impingement of the NACA 651-208 air- 
foil are compared with the results for the NACA 65,-cle airfoil. Under 
similar conditions of operation, the NACA 654-208 airfoil collects less 
water than the NACA 654-212 airfoil. The extent of impingement on the 
upper surface of the 651-208 airfoil is much less than on the upper sur- 
face of the 654-012 airfoil, but on the lower surface the extents of 


impingement are about the same. 


INTRODUCTION 


As part of a comprehensive research program directed toward an 
appraisal of the problem of ice prevention on high-speed aircraft, an 
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investigation of the impingement of cloud droplets on airfoils and other 
aerodynamic bodies has been undertaken at the NACA Lewis laboratory. The 
investigation includes a study of the extent of impingement on low-drag 
airfoils and the rate of droplet impingement per unit area of the airfoil 
surface affected. Previous investigators have calculated the water- 
droplet trajectories for cylinders (refs. 1 to 6) and for Joukowski air- 
foils (refs. 7 апа 8). Ап empirical method for determining area, rate, 
and distribution of water-droplet impingement on alrfoils of arbitrary 
sections is presented in reference 9. The method is more firmly estab- 
lished for 15-percent-thick airfoils resembling Joukowski airfoil sec- 
tions than for low-drag airfoils, because the basic data used in devel- 
oping the empirical method were obtained for four Joukowski airfoil sec- 
tions but for only one low-drag section. Some impingement data for an 
МАСА 651-212 airfoil, which ів a iZ-percent-thick low-drag section, are 


presented in reference 10. Recent developments in high-speed aircraft 
necessitate further water-droplet trajectory studies on Low-drag air- 
foils, particularly for thin sections, in order to determine the effect 


of thickness ratio on droplet impingement. 


The studies presented in this report are for 8-percent- and 12- 
percent-thick wings designated as NACA 654-208 and 651-212 airfoils, 


respectively, each placed at an angle of attack of 49. The results pre- 
sented are applicable to the NACA 651-208 airfoil and the NACA 654-212 


airfoil under the following conditions: chord lengths from 2 to 20 feet; 
altitudes from 1000 to 35,000 feet; airplane speeds from 150 miles рег 
hour to the flight critical Mach number; droplet diameters from 5 to 
100 microns; and an angle of attack of 49. The flight critical Mach num- 
ber is defined as the lowest flight speed which results in sonic veloc- 
ity at some location on the airfoil. 
SYMBOLS 

The following symbols are used in this report: 
а droplet radius, ft 
Ch drag coefficient for droplets 


pressure coefficient 


droplet diameter, microns (micron = 2.28»1079 ft) 


Qi 


2 


а У 
К inertia parameter, 2 Fw 


; Gimensionless 
9 uL 


L airfoil chord length, ft і 
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free-stream Mach number 
absolute pressure, in. Hg 


local Reynolds number with respect to droplet, 2асау/ы, 
dimensionless 


free-stream Reynolds number with respect to droplet, 2араУ/ы, 
dimensionless 


distance from element of vortex sheet to point in flow field, 
ratio to chord length 


distance on surface of airfoil measured from leading-edge chord 
point, ratio to chord length 


most probable icing temperature (fig. 15), ОВ 

time, sec 

flight speed, mph 

local air velocity, ratio to free-stream velocity 
free-stream velocity, ft/sec 

local droplet velocity, ratio to free-stream velocity 


magnitude of local vector difference between velocity of droplet 
and velocity of air, ft/sec 


rate of water impingement per unit span of airfoil, lb/(hr)(ft span) 


local rate of water impingement, lb/(hr)(sq ft) 
liquid-w&ter content in cloud, g/cu m 
rectangular coordinates, ratio to chord length 
angle of attack, deg 

local impingement efficiency, dimensionless 
vortex strength, dimensionless 

ratio of specific heats, 1.4 


viscosity of air, slug/(ft) (sec) 
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Ел coordinate points on airfoil, ratio to chord length 


р density, slug/cu ft 

T time scale, tV/LS dimensionless 
Subseripts: 

8, air 

1 lower airfoil surface 
m maximum 

в airfoil surface 

u upper airfoil surface 
У vortex 

w water 

x horizontal component 
y vertical component 


free stream 


Prime superscripts refer to coordinate system coincident with geometric 
chord line 


ANALYSIS 


As an airfoil moves through a cloud, the interception of the cloud 
droplets by the airfoil is dependent on the physical configuration of the 
airfoil, the flight conditions, and the inertia of the cloud droplets. 

In order to obtain the extent of impingement and the rate per unit area 
of droplet impingement on an airfoil, the cloud-droplet trajectories with 
respect to the airfoil must be determined. The differential equations 
that describe the droplet motion in a two-dimensional flow field have 
been derived in reference 4 and are presented herein in the following 
form: | 


x В. ) 
ат 27% Ux Ух 
il 
dv све ( ) 
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where 
p acy 


а W 2 
KEŞ I (2) 


and the Reynolds number Re is obtained in terms of the free-stream 
Reynolds number so thet 


Š 2 
(е) - (uy - vy)? (uy - vy)? (8) 


The differential equations (1) state that the motion of a droplet 
is governed by the drag forces imposed on the droplet by the relative 
motion between the droplet and the air moving along the streamlines around 
the airfoil. The droplet momentum tends to keep the droplet moving in 
a straight path, while the drag forces tend to force the droplet to 
follow the streamlines. For very small droplets and low speeds, the drag 
forces are much greater than the inertia forces and little departure from 
the streamlines occurs; whereas, for large droplets and high speeds, the 
inertia forces tend to overcome the drag forces and the droplets depart 
considerably from the streamlines and follow a path more nearly in the 
direction established by the free~stream velocity. Іп accordance with 
equations (1) and the definition of the parameter K in equation (2), 
for a given size and configuration of airfoil the trajectories depend on 
the radius of the droplets, the airspeed, and the air viscosity as first- 
order variables. The trajectories also depend on the geometry of the 
airfoil involved and its angle of attack, in that these two variables 
determine the magnitude of the component velocities of the air u, and 
Чу everywhere in the flow field. 


The component air velocities were determined by a vortex substi- 
tution method that requires a knowledge of the pressure distribution on 
the surface of the airfoil. The pressure distribution was obtained from 
wind-tunnel data taken at the NACA Ames laboratory. The fundamentals of 
the method for calculating the Local perturbation velocities in a two- 
dimensional incompressible flow field ahead of an airfoil are established 
in reference ll. The procedure for the calculations is presented in 
appendix A. The computations were performed with electronic calculating 
machines employing punched cards. Although the compressibility of air 
at high speeds is not accounted for with the vortex substitution method, 
results presented in reference 5 show that the effect of the compressi- 
bility of air on the trajectories of droplets is negligible up to the 
critical Mach number of the airfoil. Because of these studies, the 
results presented herein are considered applicable up to the flight 
critical Mach number. 
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Assumptions that have been necessary in order to solve the prob- 
lem are; 


(1) At a large distance ahead of the airfoil (free-stream condi- 
tions) the droplets do not move with respect to the air. 


(2) No gravitational force acts on the droplets. 
(3) The droplets are always spherical and do not change in size. 


The first two assumptions are valid for droplets smaller than drizzle or 
rain drops, because the inertial forces of the drops are much greater 
than the gravitational force. The assumptions are usually also valid for 
failing rain drops, since the airplane velocity is usually much greater 
than the drop velocity caused by gravitational force. Preliminary calcu- 
lations show that the third assumption is valid for the order of accuracy 
usually required in the design of equipment for the protection of air- 
erat Gs 


METHOD OF SOLUTION 


The differential equations of motion (eq. (1)) are difficult to 
solve by ordinary means, because the values of the velocity components of 
the air and the term containing the coefficient of drag are not known 
until the trajectory is traced. These values are determined as the tra- 
jectory of a droplet is developed, since the magnitudes depend on the 
position of the droplet in the flow field. Simultaneous solutions for 
the two equations were obtained with a mechanical analog constructed at 
the Lewis laboratory for this purpose. The answers were obtained in the 
form of plots of the droplet trajectories with respect to the airfoil. 
The coefficient of drag Cp for the droplets, required in equations (1), 


was obtained from tables in reference 3. 


The equations of motion (eq. (1)) were solved for the following five 
values of the parameter K: 1/100, 1/50, 1/10, 1/5, and 1. For each 
value of the parameter K, a series of trajectories was computed for each 
of three values of free-stream Reynolds number Reg: 16, 256, and 1024 


(where necessary, Бес = 64 was also computed). (A graphical procedure 


for translating the dimensionless parameters used in this report into 
terms of airplane speed, chord length, altitude, and droplet size is pre- 
sented in appendix B.) Each series of trajectories encompassed the air- 
foil with a trajectory that was tangent to the upper surface of the air- 
foil and with а trajectory that was tangent to the lower surface of the 


airfoil. The upper and lower tangent trajectories started at free-stream | 


conditions at distances Yo,u and М respectively, from the x-axis 
of the rectangular coordinate system (fig. 1). The geometric chord line 
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of the airfoil is oriented at an angle of 4° with the x-axis of the rec- 
tangular coordinate system, and the leading edge is placed at the origin 
of the coordinates. At an infinite distance ahead of the airfoil, the 
uniform air flow carrying the cloud droplets is assumed to be approaching 
the airfoil from the negative x-direction, and parallel to the x-axis. 
All distances are dimensionless, because they are ratios to the airfoil 
chord length L, which is the unit of distance. 


Before the integration of the equations of motion could be performed 
with the analog, the initial velocity of the droplets had to be deter- 
mined. As postulated in the assumptions, at an infinite distance ahead 
of the airfoll all the droplets have vertical and horizontal components 
of velocity relative to the airfoil that are the same as those of the 
free-stream air. At finite distances ahead of the leading edge of the 
airfoil, the droplets have velocity components varying between those per- 
taining to the free stream at infinity and those pertaining to the local 
streamlines. At 5 chord lengths ahead of the airfoil, the air stream- 
lines deviated from free-stream conditions by less than the expected 
accuracy of the analog; therefore, this point was assumed to be suffi- 
ciently ahead of the airfoil so that the velocity components of the air 
streamlines could be assigned to the droplets. Since the time required 
to trace each trajectory from 5 chord lengths ahead of the airfoil to the 
airfoil surface was prohibitive, plotting by the analog was started at 
l chord length ahead of the airfoil leading edge. The starting condi- 
tions at this point were determined by calculating a sample trajectory 
that started st x = -5 for each of the five values of К studied. А 
preliminary study showed that while a droplet is approaching a position 
1 chord Length ahead of the airfoil, the amount of deviation of the 
droplet trajectory from the air streamline on which the droplet had 
started depends only on the value of K and not on the starting value of 
у at х= -5, nor on tbe value of Reo, provided the values of y and 


Re, are within the region of interest in the problem. The sample tra- 


Jectories for each value of К studied were calculated from x = -5 to 
х = -l in order to determine the starting values of droplet velocity and 
y-ordinate at x = -1; thus, the final results were the same as if each 
trajectory were calculated from 5 chord lengths ahead of the airfoil 
leading edge. 


RESULTS AND DISCUSSION 


The series of trajectories computed for each combination of values 
of K and Бе, studied permits evaluation of the area, the rate, and 


the distribution of water-droplet impingement on the NACA 654-208 and the 
NACA 65,-212 airfoil sections at an angle of attack of 4°. The amount of 
water collected and the area, or extent, of impingement are determined by 
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wt 


the limiting trajectories, which are tangent to the airfoil. All droplets 
having trajectories between the tangent trajectories will strike the air- 
foil, whereas all other droplets will miss the airfoil. Тһе tangent ы 
trajectories determine the maximum, or total, rate of droplet impinge- 
ment, because the amount of water-droplet impingement on the whole wing 
is governed by the spacing between the tangent trajectories (Yo u ` Yo 7? 

2 2 


fig. 1) at a large distance ahead of the airfoil. Тһе manner in which 
all the droplets collected on the surface are distributed over the area 
of impingement is determined by the behavior of the intermediate tra- 
jectories that are bounded by the tangent trajectories. 


The results are often presented herein as functions of the parameter 
K, which has been called the inertia parameter because its magnitude 
directly reflects the external force required on a droplet to cause a 
deviation from its original line of motion. For large values of K 
(1.е., K>1), which correspond, for example, to droplets larger than 
50 microns in diameter moving toward a 5-foot-chord airfoil at 400 miles 
per hour, the droplet trajectories deviate by only small amounts from 
straight lines. For values of K less than 1/50, which correspond, for 
example, to droplets less than 12 microns in diameter moving at less than 
500 miles per hour toward an airfoil section with а 12-foot chord, the 
droplet trajectories more nearly coincide with the air streamlines. 


Rate of Water Interception ü 


In flight, the rate of total water interception (1b/(hr)(ft span)) 
is determined by both the tangent droplet trajectories and the liquid- 
water content in the cloud. The airfoil speed, the size of the airfoil, 
and the droplet size are the principal variables affecting the tangent 
droplet trajectories. The spacing between the upper and lower tangent 
trajectories, which is a measure of the rate of total water interception, 
is given in figure 2 for both the 651-208 and the 651-212 airfoils. 


The rate of total water interception per unit span of airfoil on that 
portion of the airfoil surface bounded by the upper and the lower tangent 
trajectories (fig. 1) can be calculated from the information in figure 2 


and the following relation: 


Wn = 0.529 (уо ц - Yo, 1) LUw (5) 


where the flight speed U is in miles per hour and the liquid-water 
content w is in grams per cubic meter. 


For both airfoils, the rate of water interception is decreased as 
YE is increased, particularly for values of 1/K larger than 1. Ап 
increase in Reg decreases slightly the rate of water interception. 


The following example illustrates the use of figure 2 and equation (5) 
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in comparing the two airfoils in terms of conventional flight units. The 
values of 1/К ж 21 and Reg = 171 represent the conditions for airfoils 


with lZ-foot chords traveling at 400 miles per hour through a cloud com- 
posed of droplets 17 microns in diameter at an altitude of 10,000 feet 
(1/K and Вер based on most probable icing temperature, appendix B). 


The corresponding value for (Yo,u - yo,1) required for equation (5) is 
0.017 for а 651-208 airfoil as compared with 0.021 for а 654-212 airfoil. 


If the same airfoils are considered in flight at 200 miles per hour 
(droplet size, liquid-water content, and altitude not varied), the values 
of 1/K and Reg change to 4@ and 85, respectively, and the value of 


(You - У0,1) changes to 0.015 for the 8-percent-thick airfoil and to 


0.015 for the lZ-percent-thleck airfoil, a decrease of 24 and 28 percent, 
respectively, from the corresponding values associated with flight at 

400 miles per hour. The decrease in rate of water impingement (eq. (5)) 
is 62 percent for the 8-percent-thick airfoil and 64 percent for the 
le-percent-thick airfoil. The effect of speed оп the rate of water 
impingement is large, because the spacing between the two tangent trajec- 
tories increases with the increase in flight speed and the speed also 
appears directly in equation (5). 


The rate of water interception is less on the 8-percent-thick airfoil 
than on the lZ-percent-thick airfoil over the complete range of free- 
stream Reynolds number and inertia parameter studied (flight speed, alti- 
tude, droplet sizes encountered in clouds, and airfoil size) except for 
combinations of 1/K>65 and Re j<32, which correspond, for example, to 


flight with an airfoil chord length larger than 5 feet at a speed of less 
than 150 miles per hour through clouds at 15,000 feet altitude composed 
of droplets less than lO microns in diameter. 


The variation of rate of water interception with airfoil speed is 
summarized for an altitude of 20,000 feet in figure 3, in which the 
ordinate Wu is the total rate of water impingement per foot span of 


airfoil per unit liquid-water content (g/eu m) in the cloud. (Fig. 3 
is based on most probable icing temperature, appendix B). The total rate 
of water impingement can be obtained as a product of the results in fig- 
ure 5 and the liquid-water content existing in the cloud. The values in 
figure 3 are for flight through a uniform cloud composed of droplets 15, 
20, 30, and 40 microns in diameter, with several chord lengths ranging in 
value from 2 to 20 feet considered. 


The effect of altitude on the rate of water impingement is shown in 
figure 4 for droplets smaller than 40 microns in diameter. The values 
for rate of water impingement given in figure 3 can be corrected to a 
first order of approximation for the altitude effect with the use of fig- 
ure 4. For example, the rates of water impingement given in figure 3(c) 
can be corrected for use at 10,000 feet altitude by reducing the values 
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of rate of water impingement given in figure 3(с) by 7 percent, because 
for 50-шісгоп droplets the rate of water impingement at 10,000 feet is 
approximately 7 percent less than at 20,000 feet (fig. 4). A change in 
altitude of 10,000 feet for any given droplet size shown in figure 4 does 
not produce as large an effect on the rate of water impingement as does 
а 15-регсепі change in droplet size at a given altitude. Тһе effect of 
moderate changes in altitude may be ignored, because the droplet size 
and the liquid-water content of clouds are seldom known with sufficient 
accuracy to permit the rate of water collection to be calculated within 
10 percent. Therefore, within the practical limits of application, the 
results of figure 3 can be used over a wide range of altitudes. 


The effect of wing taper on the rate of water impingement, as can 
be obtained from the results in figure 3, is valid, provided that for 
each section of span considered the taper is small enough to permit the. 
approximation of two-dimensional flow over the section. For a 651-208 


airfoil section with a 14-Ғоо% chord at the root section and a 4-foot 
chord at the tip section moving at 400 miles an hour at an altitude of 
20,000 feet through a cloud composed of droplets 20 microns in diameter 
and containing a liquid-water content of 1 gram per cubic meter, the 

rate of total water impingement is 35 pounds per hour per foot span at 

the root section and 17 pounds per hour per foot span at the tip section. 
A lZ-percent-thick airfoil of the same size and section subjected to the 
same flight and atmospheric conditions as given for the 8-percent-thick 
airfoil collects approximately 28 percent more water at both the root and | 
tip sections. 


Extent of Impingement 


A knowledge of the chordwise extent of impingement on the upper and 
lower wing surfaces is necessary for the design of anti-icing equipment. 
The limit of impingement is determined by the point of tangency on the 
airfoil of the tangent trajectories. The rearward limits of impingement 
on the upper surface are shown in figure Ба), and those on the lower 
surface, in figure 5(b). Тһе distances S, and 8) are measured оп the 


surface from the point of intersection of the geometric chord line with 
the leading edge (fig. 1) in terms of the chord length. 


Ihe limit of impingement on the 651-208 airfoil is compared in fig- 
ure 5 with that on a 65,-212 airfoil. The impingement extends much farther 


back over the top surface of the 12-percent-thick airfoil than over the 
top surface of the 8-percent-thick airfoil (fig. 5(a)). The impingement 
extends farther back along the bottom surface of the thicker airfoil over 
most flight and atmospheric conditions, but is farther back on the lower 
surface of the 8-percent-thick airfoil for combinations of flight and 
atmospheric conditions involving 1/K values smaller than 1.5 and larger 
tham 50. 
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The rearward limit of impingement along the upper and lower sur- 
faces, respectively, is summarized in figures 6 and 7 for the same 
speeds, chord lengths, droplet sizes, and altitude ав given in figure 5.: 
The extent of impingement on both the upper and lower surfaces increases 
with increasing speed and with decreasing chord length. It is much 
greater along the lower surface than along the upper surface. As an 
example, for the 65,-208 airfoil with a 14-foot chord traveling at 


300 miles per hour through a cloud composed of droplets ZO microns in 
diameter, the extent of impingement on the upper surface is 0.006 chord 
(fig. 6(b)), and on the lower surface is 0.071 chord (fig. 7(b)). Again, 
as shown in figure 5, the extent of impingement on the 8-percent-thick 
airfoil is less than on the lZ-percent-thick airfoil, except on the lower 
surface for combinations of low flight speed and large chord length and, 
particularly, with the smaller droplets considered in figure 7(8). These 
combinations of low flight speeds and large chord lengths given in fig- 
ure 7(a) correspond to conditions involving values of 1/K larger than 
50 shown in figure 5(b). The combination of conditions involving values 
of 1/К<1.5, for which the extent of impingement on the lower surface of 
the 8-percent-thick airfoil is also greater than on the lower surface of 
the 12-percent-thick airfoil, is given in figure 7(c) and 7(d). Values 
of 1/K<1.5 are involved, for example, with combinations of airfoil 
chord lengths smaller than 4 feet and flight speeds greater than 400 miles 
per hour for a droplet size of 40 microns, as shown in figure 7(d). 


impingement Distribution Along Airfoil Surface 


In the design of thermal anti-icing systems based either on the prin- 
ciple of maintaining the water in the liquid state or of completely 
evaporating the impinging water, a knowledge of the distribution of the 
water impingement along the airfoil surface is required. The distribu- 
tion сап be obtained if the starting ordinate yọ at infinity of а drop- 


let trajectory is known with respect to the point of impingement on the 
surface. The starting and ending positions of the trajectories are 
shown in figure 8 for the three values of free-stream Reynolds number 
studied. For each value of Reg, curves for several values of 1/K are 


given. 


The amount of water impinging between any two given points on the 
airfoil surface may be found by applying the results given in figure 8 
in the relation 


М = о. 52909 (у, 5 - Yo,2) (8) 


For example, the amount of water impinging between the -0.05-chord point 
and the -O.10-ehord point on the surface of a 12.5-foot-chord airfoil 
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moving at 400 miles per hour at ап altitude of 10,000 feet through a 
cloud composed of droplets 25 microns in diameter (1/К = 10, Reg Z 256) 


is determined by obtaining the values of Уо,2 at the -O.lO-chord point 
land yo at the -O.O5-echord point. The values of У0,1 and 30:5 
required in equation (8) are -0.074 and -0.078 for the 8-percent-thick 


airfoil and -0.071 and -0.076 for the lZ-percent-thick airfoil, respec- 
tively (fig. 8). 


The values of (У, LI 1) obtained from the end points of each 
2 2 


curve in figure 8, are the same as the values given in figure 2. Те 
values of усу for 1/K = О (not shown in fig. 8) are -0.0644 at 


8 = -0.654 chord for the 8-percent-thick airfoil and -0.0825 at 
5 = -0.538 chord for the lZ-percent-thick airfoil. These values Гог 
1/K = О apply for all values of free-stream Reynolds number Reg. 


Local Rate of Droplet Impingement 


The local rate of droplet impingement per unit area of airfoil sur- 
Pace can be determined from the expression 


= = 0.529UvB (9) 


Ше = 0.329Uw 
This equation is related to equation (8), with proper considerations for 
the fact. that yo and 5 are based on the wing chord L. The values 


of the local impingement efficiency В as a function of the airfoil dis- 
tance S are given in figures 9 and 10. These values were obtained from 
the slopes of the curves in figure 8. The values of 8 for the NACA 
654-208 airfoil are given in figure 9, and for the 651-212 airfoil, in 


figure lO. For both airfoils, the maximum rate of local impingement 
occurs between the air stagnation line and the leading edge. The air 


stagnation line is located at 5, = 0.013 on the 65, -208 airfoil and at 


8) = 0.008 оп the 65,-212 airfoil (at 4° angle of attack). 


Because of the geometry of the airfoil and the manner in which the 
droplets approach the airfoils in the neighborhood of the stagnation line 


and the leading-edge line, the curves in Pigure 8 are not all well defined 


between 5 = -0.01 and 5 = O. Since the values of В are obtained from 
the slopes of the curves in figure 8, the maximum values of В in fig- 
ures 9 and 10 also are subject to some question. The possible error in 
the maximum value of В is estimated to be +20 percent for the curves in 
figure 9 and +15 percent for the curves in figure 10. The difference in 
the possible error is due to the difference in the bluntness of the two 
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airfoils. This possible error is not considered very serious, because 
only a small portion of the total water impinging on the airfoil requires 
redistribution if the maximum value of В is changed by аз much as the 
maximum possible error. The total area under the curves should not be 
changed when а change in the maximum value of В is made, because the 
total amount of water impinging, determined by oes - Yo 1) and by 


S 
ДЕ 8 48, is independent of the manner in which the water is distributed 
near the leading edge. Апу judicious refairing of the curves in Тір- 

ures 9 and 10 between 5 ж -О.О1 and 5 = 0, which is the region of 
doubt, to &ecount for as mueh as the largest possible error in the maxi- 
mum value of В will rearrange only a small area under the present 

curves (total area must not be changed) and thus affect the redistribu- 
tion of only а small amount of the total water. 


The maximum local rate of impingement shown for these airfoils in 
figures 9 and 10 coincides very nearly with the location of the parting 
strip recommended in references 12 and 13 for use in cyclical thermal 
de-icing systems. 


CONCLUDING REMARKS 


The data presented in figures 2 to 10 apply directly to flights in 
clouds composed of droplets that are all uniform'in size. The water 
droplets in & cloud, however, are not necessarily uniform in size; the 
extent of impingement is always determined by the largest droplets pres- 
ent in sufficient number to represent a significant portion of the total 
water present in the cloud. The local rates of droplet impingement are 
also determined by the droplet-size distribution patterns present in the 
cloud. For flights in clouds composed of droplets that are not uniform 
in size, the В curves of figures 9 and 10 must be altered to conform 
with the weighted basis of the droplet-size distribution. A detailed 
procedure for weighting the impingement of droplets existing in nonuni- 
form clouds is presented in reference 4. 


The data presented in this report apply directly to nonswept wings 
of high aspect ratio. As previously stated, the data also apply to wings 
with small amounts of taper. А method for extending the impingement 
calculations for nonswept wings to swept wings is presented in refer- 
ence 14. The method presented in reference 14 permits the application 
of the data presented herein to a wing with any angle of yaw provided the 
airfoil section normal to the leading edge is an NACA 65,-208 or ап NACA 
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651-212 airfoil. ТЕ the airfoil section on the aircraft under considera- 
tion is laid out parallel to the free-stream direction, a transformation 


of the section is necessary, as described in reference i4. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 


Cleveland, Ohio, January 30, 1955 
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APPENDIX A 


METHOD USED TO CALCULATE INCOMPRESSIBLE FLOW FIELD AROUND AIRFOIL 


The velocities in the two-dimensional flow field were calculated by 
distributing a sheet of vortices on the airfoil surface of such strength 
that the velocities on the surface caused by the vortices were the same 
as the velocities measured in a wind tumnel. The principles upon which 
this method is based are established in reference 11. 


Ihe velocity at tbe surface of an airfoil can be determined from a 
knowledge of the pressure coefficient Ср and the free-stream Mach 


number M with the aid of the following expression: 


(A1) 


Ihe pressure coefficients for & large number of points on the surface of 
the NACA 651-208 and 651-212 airfoils at an angle of attack of 4° were 
obtained from wind-tunnel data taken at the NACA Ames laboratory for 
several free-stream Mach numbers. The surface velocities used to cal- 
culate the flow field are shown in figure 11 for а Mach number of 0.2. 
Тһе flow fields at other Mach numbers were not calculated, because the 
results presented in reference 5 show that the effect of the compressi- 
bility of the air on the droplet trajectories is negligible. 


Тре velocity at а point in a flow field caused by an element of the 
vortex sheet of strength ЛГ placed a distance г (fig. 12) from the 
point is 


If an element of vortex sheet of strength 


is placed on an increment AS of the airfoil at the 190 position оп the 
airfoil surface, the velocity caused only by the ith section of the 
airfoil is 


(ug^8) з 


V алг+ 
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at a point in the flow field at a distance x, from the ¿th section 


(fig. 12). The local components of the perturbation velocity, at a 
point in the flow field, caused by 300 vortex elements distributed on 
both the upper and lower surfaces of the airfoil are 


500 


OO us Ё-х' 
ум” "E A 28 


The coordinate system shown in figure 12 differs from that shown in fig- 
ure 1, in that in figure 12 one of the coordinates coincides with the 
geometric chord line. The direction of the free-stream velocity is 
toward the airfoil at an angle a with respect to -x'. The horizontal 
and vertical components of the local velocity u,: and Uyt, respec- 
tively, are obtained by adding V cosa to Их! iy and V sina to 
uy',v: The coordinate system shown in figure 12 was used because of its 
adaptability to available calculating equipment and procedures. The 
trajectories were solved in the primed coordinate system shown in fig- 
ure 12 and later transformed graphically to the system presented in fig- 
ure 1. 


A total of 300 vortex elements were used on the airfoil with a much 
denser distribution on the forward section than beyond the 50-percent- 
chord point. Equations (A2) were solved with electronic calculating 
machines. Approximately 300 points were computed in the flow field out 
to L chord length ahead of the airfoil in the region of interest with 


regard to computing the trajectories of droplets that strike the airfoil. 


Between 1 and 5 chord lengths ahead of the airfoil leading edge, the 
flow-field velocity components were approximated by assuming that the 
flow was caused by a single vortex located on the airfoil chord Line 
25-percent chord inward from the leading edge. The strength of this 
vortex was determined by the requirement that at хі -l the vertical 
velocity caused by this single vortex must be the same as that computed 


with equation (A2). 
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APPENDIX B 


GRAPHICAL PROCEDURE FOR TRANSLATION OF PRACTICAL FLIGHT 
CONDITIONS IN TERMS OF DIMENSIONLESS PARAMETERS 


A graphical procedure is presented to ада in the translation of 
airplane speed, chord size, altitude, and droplet diameter into terms of 
the dimensionless parameters K and Rep used in this report. A 
solution of equation (2) is presented in figure 13 for two altitudes. 
For given droplet diameters in- microns and ratios of the chord length 
in feet to the flight speed in miles per hour, the reciprocal of the 
inertia parameter can be determined at altitudes of either 10,000 or 
50,000 feet from figure 15. Altitude does not appreciably affect the 
value of 1/K, as can be noted from a comparison of values in fig- 
ure 13(a) with those in figure 13(b). 


An airfoil with а 12-foot chord length at а flight speed of 
400 miles per hour and an altitude of 10,000 feet passing through a | 
cloud composed of droplets all of which are 17 microns in diameter will' 
be used as an example in the graphical procedure to interpret practical 
flight units into terms of the dimensionless parameters. The value of 
l/K is obtained from figure 15(а) for the value of = = = = 0.0500 
and droplet diameter d = 17. The value of 1/K obtained from fig- 
ure 13(а) is 21. 


The free-stream Reynolds number for different altitudes may be 
obtained from figure 14. The product of the droplet diameter in microns 
and the flight speed in miles per hour must be known. The Reynolds 
number is a function of the air density, which depends on the pressure 
and the temperature at the altitudes considered. The pressure used to 
calculate the air density was taken from tables of NACA standard atmos- 
pheric pressure at various altitudes, but the temperature was based on 
the most probable icing temperature at various altitudes. The most 
probable icing temperature was obtained from approximately 300 icing 
observations (ref. 15) and is presented in figure 15. For the example 
under consideration, the product of the droplet diameter and the flight 
speed is (17)(400) = 6800. The value of Rep, obtained from figure 14, 
is 171. 


The following relations are presented for use when the degree of 
accuracy required is not attainable with the graphical procedure. The 
values for the viscosity ы should be obtained from figure 16; these 
values are based on the most probable icing temperature of figure 15. 
The charts of figures 13 and 14 are based on the most probable icing 
temperature and viscosity. 
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4.815 X 10 


7.662 X 10? E 


р 
ра = 0.0412 т. 
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where the units are given in the llst of symbols. (The density of water 
was assumed to be 62.46 lb/cu ft and the acceleration due to gravity, 
52.17 ft/sec.) 


REFERENCES 


l. Glauert, Muriel: A Method of Constructing the Paths of Raindrops of 
Different Diameters Moving in the Neighborhood of (1) & Circular 
Cylinder, (2) an Aerofoil, Placed in a Uniform Stream of Air; and 
a Determination of the Rate of Deposit of the Drops on the Surface 
and the Percentage of Drops Caught. Е. & M. No. 2025, British 
A.R.C., 1940. 


2. Ranz, W. E.: The Impaction of Aerosol Particles on Cylindrical and 
Spherical Collectors. Tech. Rep. No. 3, Eng. Exp. Station, Univ. 
Ill., March 31, 1951. (Contract No. AT(30-2)-28, U.S. Atomic 
Energy Commission.) 


3. Langmuir, Irving, and Blodgett, Katherine B.: A Mathematical Inves- 
tigation of Water Droplet Trajectories. Tech. Rep. No. 5418, Air 
Materiel Command, AAF, Feb. 19, 1946. (Contract No. W-33-038-ac- 
9151 with General Electric Co.) 


4. Brun, Rinaldo J., and Mergler, Harry W.: Impingement of Water Drop- 
Lets on a Cylinder in an Incompressible Flow Field and Evaluation 
of Rotating Multicylinder Method for Measurement of Droplet-Size 
Distribution, Volume-Medi&n Droplet Size, and Liquid-Water Content 
іп Clouds. ЛАСА TN 2904, 1953. | 


5. Brun, Rinaldo J., Serafini, John S., &nd Gallagher, Helen M.: | 
Impingement of Cloud Droplets on Aerodynamic Bodies as Affected by 
the Compressibility of Air Flow Around the Body. NACA TN 2903, 
1953. 


с942 


3912 


CN-3 back 


NACA TN 2952 19 


6. Brun, Edmond, Caron, Robert, et Marcel, Vasseur: Introduction Á 
l'étude de la Mécanique des Suspensions. G.R.A. Rapport Tech. 
No. 15, Recherches Aéronautiques (Paris), 1945. 


1. Bergrun, Norman R.: A Method for Numerically Calculating the Area 
and Distribution of Water Impingement on the Leading Edge of an 
Airfoil in a Cloud. МАСА TN 1597, 1947. 


8. Guibert, A. G., Janssen, E., and Robbins, W. M.: Determination of 
Rate, Area, and Distribution of Impingement of Waterdrops on 
Various Airfoils from Trajectories Obtained on the Differential 
Analyzer. КАСА RM 9А05, 1949. 


9. Bergrun, Norman R.: An Emperical Method Permitting Rapid Determina- 
tion of the Area, Rate, and Distribution of Water-Drop Impingement 
on an Airfoil of Arbitrary Section at Subsonic Speeds. NACA TN 
2476, 1951. 


lO. Brun, Rinaldo J., Serafini, John 5., and Moshos, George 4.: 
Impingement of Water Droplets on ап NACA 651-212 Airfoil at an 
Angle of Attack of 49. МАСА RM E52B12, 1952. 


11. von Mises, Richard: Theory of Flight.  MeGraw-Hill Book Co., Inc., 
First ed., 1945. 


l2. Gray, V. HE., Bowden, D. T., and von Glahn, U.: Preliminary Results 
of Cyclical De-icing of a Gas-Heated Airfoil. NACA RM Е51429, 
1952. 


13. Lewis, James P., and Bowden, Dean T.: Preliminary Investigation of 
Cyclic De-Icing of an Airfoil Using an External Electric Heater. 
NACA RM E51J30, 1952. 


14. Dorsch, Robert G., and Brun, Rinaldo J.: A Method for Determining 
Cloud-Droplet Impingement on Swept Wings. МАСА TN 2931, 1953. 


15. Hacker, Paul T., and Dorsch, Robert G.: А Summary of Meteorological 
Conditions Associated with Aircraft Icing and a Proposed Method of 
Selecting Design Criterions for Ice-Prevention Equipment. МАСА 
UN 20593. 1902» 


— f 5 (опері л chord line 


pper tangent trajectory 40 7-- u 


Lower tangent trajectory 


X в - сь 


Figure 1. ~ Droplet trajectories with respect to airfoil. 
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Figure 2. - Spacing between upper and lower tengent trajectories at fres-streanm conditions, Angle of attaok, 4°. 
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Figure 3. - Rate of water impingement. Angle of attack, 4°; alti- 
tude, 20,000 feet; most probable icing temperature, -119 Е. 
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Figure 5. - Continued. Rate of water impingement. Angle cf 
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Figure 8. - Trajectory starting ordinates as function of point of impingement on airfoil вигТасе. 
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Figure 9. - Continued. 
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(Based on most probable 


Figure 16. - Air viscosity as function of atmospheric altitude, 
icing temperature.) 


